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Abstract —We model the self-interference in a multiple input 
single output (MISO) full-duplex two-way channel and evaluate 
the achievable rate region. We formulate the boundary of the 
achievable rate region termed as the Pareto boundary by a 
family of coupled, non-convex optimization problems. Our main 
contribution is decoupling and reformulating the original non- 
convex optimization problems to a family of convex semidefi- 
nite programming problems. For a MISO full-duplex two-way 
channel, we prove that beamforming is an optimal transmission 
strategy which can achieve any point on the Pareto boundary. 
Furthermore, we present a closed-form expression for the optimal 
beamforming weights. In our numerical examples we quantify 
gains in the achievable rates of the proposed beamforming over 
the zero-forcing beamforming. 

I. Introduction 

A node in full-duplex mode can simultaneously transmit and 
receive in the same frequency band. Therefore, the wireless 
channel between two full-duplex nodes can be bidirectional, 
having the potential to double the spectral efficiency when 
compared to the half-duplex network. Due to the proximity 
of the transmitters and receivers on a node, the overwhelm¬ 
ing self-interference becomes the fundamental challenge in 
implementing a full-duplex network. The mitigation of the 
self-interference signal can be managed at each step of the 
communication network by passive and active cancellation 
methods (T) ■ In recent results the feasibility of the sin¬ 

gle input single output (SISO) full-duplex communication has 
been experimentally demonstrated. However, the performance 
is limited by the residual self-interference which is considered 
in 0>0-@ to be induced by the imperfection of the transmit 
front-end chain. 

The bottleneck from imperfect transmit front-end chain has 
motivated recent research in full-duplex channel with transmit 
front-end noise. The performance of the SISO full-duplex two- 
way channel has been thoroughly analyzed in 0 > 0 - The 
multiple input multiple output (MIMO) full-duplex two-way 
channel with transmit front-end noise is considered in 0, @ 
(in 0, |6| termed as MIMO full-duplex bidirectional channel). 
In |5| Vehkapera et al. studied the effect of time-domain 
cancellation and spatial-domain suppression on the channel, 
while in |6| Day et al. derived the lower bound of achievable 
sum-rate for the channel and proposed a numerical search for 
optimal signaling. In this paper, we focus on the multiple 
input single output (MISO) full-duplex two-way channel in 
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Fig. 1. The MISO point-to-point full-duplex network under study. The 
solid line denotes the desired channel and the dashed line denotes the self¬ 
interference channel. 


presence of the transmit front-end noise. Compared with 0, 
we derive the tight boundary of the achievable rate region for 
the channel and present the analytical closed-form solution 
of the optimal signaling. Note that the achievable rate region 
includes the achievable sum-rate as a point and provides the 
additional asymmetric performance metric. 

In this paper, we consider the optimal signaling structure 
for the MISO full-duplex two-way channel, by which all rate 
pairs on the boundary of the achievable rate region can be 
achieved. We introduce the channel model that includes trans¬ 
mit front-end noise. We leverage our model to characterize 
the achievable rate region for the full-duplex channel. The 
boundary of the region is described by a family of non- 
convex optimization problems. Rendering the computation 
tractable, we decouple the original non-convex problems to 
the family of convex optimization problems. The decoupling 
method was first developed in the field of game theory @ 
and recently introduced to communications in 0-G3- By 
employing the semi-definite programing (SDP) reformulation, 
we numerically solve the optimal signaling and prove the 
optimality of transmit beamforming. That is to say, for a MISO 
full-duplex two-way channel, all the points on the boundary 
of the achievable rate region can be achieved by restricting 
to transmit beamforming scheme. Furthermore, we derive the 
closed-form optimal beamforming weights. Finally, through 
simulations we show the achievable rate regions for the MISO 
full-duplex two-way channels and evaluate the performance 
of the traditional zero-forcing beamforming with our optimal 
beamforming. 

Notation: We use (-)t to denote conjugate transpose. For a 
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Fig. 2. The model of the MISO full-duplex two-way channel under study. 

scalar a, we use |a| to denote the absolute value of a. For a 
vector a £ C Mxl , we use ||a|| to denote the norm, a ^ to 
denote the k th element of a, Diag(a) to denote the square 
diagonal matrix with the elements of vector a on the main 
diagonal. For a matrix A £ c MxM , we use A _1 , tr(A) 
and rank (A) to denote the inverse, the trace and the rank 
of A, respectively. We use diag(A) to denote the diagonal 
matrix with the same diagonal elements as A. A >; 0 means 
that A is a positive semidefinite Hermitian matrix. We denote 
expectation, variance and covariance by E{-}, Var{-} and 
Cov{-}, respectively. Finally, C and H denotes the complex 
field and the Hermitian symmetric space, respectively. 

II. Channel Model 


/3diag (Qj) where i3 is a constant depending on the distortion 
level of the transmit front-end chain |I |. 

At receiver i , the signal of interest j ^ i is received 

along with the self-interference signal b.k.s, and the transmit 
front-end noise h^ej, he,. The power level of ^ i 

is typically much lower than that of the thermal noise Ni and 
thus can be neglected JiJ. However, h ' n e, is in the power level 
close to the signal of interest and needs to be considered for 
the analysis, since the gain of the self-interference channel ha 
may be lOOdB higher than the gain of the cross-node channel 
hji |6j. In addition to the strength, transmitters and receivers 
on a same node are relatively static, resulting in long coherence 
time of self-interference channels, thus receiver i is assumed to 
have the perfect knowledge of its own self-interference channel 
ha (5j. Note that receiver i also knows its own transmitted 
signal ,s,. Then we can eliminate the self-interference h] ? s, 
before decoding. The signal after cancellation is given by 

Yi = h) iSj + h\ iei + Ni (2) 

where h' u ei represents the residual self-interference. By defin¬ 
ing the aggregate noise term V = h'^ei + Ni, the received 
signal model can be further simplified as 

Yi = h) iSj + Vi (3) 

where Vi is Gaussian noise with zero mean and variance 
Var{VJ = ct 2 n + ph\i&mg(Qi)hii. 


We present the channel model for a MISO full-duplex 
network with two nodes as illustrated in Fig. |T| Assume 
two nodes indexed by i,j £ {1,2} share the same single 
frequency band for transmission. Each node is equipped with 
M transmit antennas and a single receive antenna. The signal 
from transmitter i is collected as the signal of interest by 
receiver j, j ^ i, while appears at its own receiver i as the 
self-interference signal. 

Fig. ^summarizes our model of the MISO full-duplex two- 
way channel. Denote the wireless channel from transmitter i 
to receiver j by the complex vector hij £ C Mxl . The signal 
at receiver i is given by 

Yi = hji(sj + ej) + h\,i(si + e^) + Ni (1) 

where Si £ C Mxl denotes the transmit signal prior to the 
transmit front-end chain at transmitter i. An additional transmit 
front-end noise e, is propagated over the same channel as 
Si. At receiver i, the thermal noise is modeled as a complex 
Gaussian noise Ni ~ CAf(0, <J 2 N ). 

The transmit front-end noise e, is induced by the imperfect 
transmit front-end chain 0-@. More precisely e, statistically 
relates to the transmit signal Si due to the limited dynamic 
range of the transmit front-end chain 0, §■ Denote the 
covariance of Si by Q, = Cov {.s,}. Note that the diagonal 
elements of Q t represent the transmit signal power. Following 
the results in we model e, as the independent 

Gaussian vector with zero mean and covariance Cov{e,} = 


III. Boundary of Achievable Rate Region 

We present the achievable rate region for the MISO full- 
duplex two-way channel and characterize the boundary points 
of this region by a family of coupled nonconvex optimization 
problems. Next, we show that the boundary points can be 
alternatively obtained by solving a family of convex optimiza¬ 
tion problems that are the results of the transformation of the 
original nonconvex optimization problems. 


A. Achievable Rate Region and Pareto Boundary 

As shown in 0, the channel model for the wireless link 
from node 2 to node 1 is equivalent to a Gaussian channel. 
It follows from the results of (B 03 that by employing a 
Gaussian codebook at node 2, we can achieve the maximum 
rate for the channel from node 2 to node 1 


RliQnQi) = log 1 + 


^ 21 ^ 2^21 


t 2 N + /3/4idiag(Q 1 )hi 


(4) 


where {Q 1 ,Q 2 ) are the given transmit covariance matrices. 
Similarly, the maximum rate for the channel from node 1 to 
node 2 is equal to 


# 2 (Qi,Q 2 ) = log 1 + 


^' 12 ^ 1^12 


o% + /?/4 2 diag(Q 2 )/i 2 


(5) 


Any rate pair (ri,r 2 ) with n < i?i,r 2 < f? 2 , is achievable 
for the MISO full-duplex two-way channel. 



























Define the achievable rate region for the MISO full-duplex 
two-way channel to be the set of all achievable rate pairs under 
the transmit power constraint Pg. 

f (n,r 2 ) : 1 

n = u { Q<ri<Ri(Qi,Q 2 )\ ■ (6) 

lo<r 2 <f? 2 (Q l5 Q 2 )J 

A rate pair is on the boundary of the rate region TZ if and only 
if it is Pareto optimal which is defined as follows (A similar 
definition can be found in 03, 03- ©)■ 

Definition 1: (Pareto optimality) A rate pair (R\,R 2 ) £ 
TZ is Pareto optimal if there does not exist another rate pair 
(f?i, f? 2 ) £ TZ such that (Ri,R 2 ) > (R*,R 2 ) and (i?i, R 2 ) 
(Rl,R 2 ) where the inequality is component-wise. 

Accordingly, the set of all boundary points of the achievable 
rate region TZ is called Pareto boundary and defined as 

TZ* = {all the Pareto optimal rate pairs (. R 2 ) }. (7) 

It is shown in © that TZ* can be derived by solving a 
family of nonconvex optimization problems: 

max /j,iRi{Q 1 ,Q 2 ) + V 2 R 2 (Qi,Q 2 ) 

(o) 

subject to tr(Qj) < Pi, Q z h 0, i = 1, 2 

where 0 < g\ < oc and 0 < p 2 < oo. However, the non¬ 
convexity of problem ([8]) implies that we must go through 
all possible transmit covariance matrices Q 1 and Q, to find 
the optimal solution for each (pi, p 2 ) pair. What is worse, 
the complexity of such exhaustive search exponentially in¬ 
creases with the dimensions of Q l and Q>, which renders the 
computation intractable E3- Next, we present an alternative 
approach more suitable for deriving the Pareto boundary for 
the MISO full-duplex two-way channel. 

B. Decoupled Optimization Problems 

The difficulty in deriving Pareto boundary TZ* is caused by 
the non-convexity and the coupled high-dimensional nature of 
problem 0. To reduce the computationally complexity, we 
need to decouple problem 0 in terms of lower-dimensional 
variables. Using the decoupling procedure in |9j, flO) , 1 13) we 
introduce an auxiliary variable z-i for node i which denotes the 
power of the signal of interest at node j i.e., Zi = h\ ] Q l h,j. 
Consider the following optimization problem for node i with 
the transmit power constraint Pg. 

min h[ l dvAg(Q i )h ll 

subject to If Q,h,, = z. t (9) 

tr(Qi) < Pi , Qi h 0 

where i, j £ {1,2} and i j. We require 


in terms of z % and F* (zfi as follows: 

' (n,r 2 ) : 


TZ 


u 


*ie[o,Pi||h 12 f 

Z2e[0,P 2 ||tl 2 l|| 2 


n = log l + 


r 2 = log 1 + 


^2 


a 2 N + pT\(zi) 
°n + /3r 2 (22) 


. 


Any rate pair in the above set TZ can be achieved by the optimal 
solution Q\ and Q* 2 for problem (|9]» with some Z\ and z 2 , thus 
TZ is a subset of the achievable rate region TZ in In Lemma 
1 we show that TZ includes the Pareto boundary TZ* for the 
region TZ. 

Lemma 1: Under the transmit power constraint Pi, any point 
on the Pareto boundary TZ* for the achievable rate region TZ 
in (|0) can be achieved by the optimal solution Q* for problem 
(|9]» with some That is to say, TZ* C TZ. 

Proof: For any point (f?{,f? 2 ) on the Pareto boundary, 
assume that it is achieved by Q\ and Q 2 . Q* is a feasible 
solution for problem (jhji with z, = z* = h[ :i Q*h l3 where 
i,j £ {1,2} and i Let i = 1, if is not an optimal 
solution for problem (j9j» i.e., /i{ 1 diag(Q{)/in > r{( 2 {) then 


R\ < log 



ofr + WM) 


— Ri, 


while 


R*2 < log 





= f? 2 . 


As (Ri,R 2 ) belongs to TZ and thus belongs to TZ, < R\ 
and R 2 < R 2 contradict to the Pareto optimality of {R\,R 2 ). 
Therefore Q\ is an optimal solution for problem 0. In the 
same way we can show that Q 2 is an optimal solution for 
problem ©• _ ■ 

We stress that the set TZ is not necessarily equivalent to the 
Pareto boundary TZ*, since TZ may include the rate pairs inside 
the region TZ. However, the relationship TZ* C TZ implies that 
any approach of obtaining the set TZ will suffice to derive the 
entire Pareto boundary TZ*. Furthermore, any result applying 
to TZ also works for TZ*. Hence, we proceed to explore the 
optimal signaling for the MISO full-duplex two-way channel 
by the study of the set TZ. 


IV. Optimal Signaling 

By solving problem 0, we show the optimality of transmit 
beamforming for the MISO full-duplex two-way channel. In 
other words, all the points on the Pareto boundary TZ* of 
the achievable rate region TZ can be achieved by transmit 
beamforming. To better understand the optimal signaling, we 
provide the closed form of the optimal beamforming weights. 


0 <Zi< max hl.Qihij = Pi\\hij \\ 2 ( 10 ) 

tr{Qi)<Pi.QihO J 

so that problem 0 always has a feasible solution. Denote the 
optimal value of problem 0 as V* ( Zi ). Then, we define a set 


A. Semidefinite Programming Reformulation 

Problem 0 is not a common optimization problem since 
the objective function includes the non-linear operator diag(-). 
By setting = h,,h[ r C = Diagd/i^| 2 ,..., |^ M) | 2 ) and 










using the equivalent relationship h\ i diag(Q i )hu = tr(C,Q, : ), 
we reformulate problem (|9]) to the semi-definite programming 
(SDP) problem as follows (See more details about SDP in 

C 3 ): 


min tr(C t Qj) 
subject to tr(AjQJ = Zi, 

tr (Qi) < Pi, Q, >: 0 


( 11 ) 


where Ci, A,- £ H M . 

The above SDP reformulation reveals the hidden convex¬ 
ity of problem 0 so that we can solve it by employing 
the well-developed interior-point algorithm within polynomial 
time. Furthermore, we can numerically characterize the Pareto 
boundary for the MISO full-duplex two-way channel in effi¬ 
ciency. 

B. Optimal Beamforming 


The optimal solutions for problem (11 1 determine the signal 


ing structure to achieve the rate pairs in the set 7Z. In Theorem 
1, we explore the rank of optimal solutions Q* for problem 


(111 where i, j £ { 1, 2} and i ^ j. 

Theorem 1: For problem GD with Pi > 0 and 0 < z t < 
Pi || hij || 2 , there always exists an optimal solution Q* with 
rank(Q*) = 1. 

Proof: See the Appendix. ■ 

Note that the transmit signal with the rank-one covariance 
matrix can be implemented by transmitter beamforming. It 
follows from Theorem 1 that all points in the set 7 Z, which 
include the entire Pareto boundary, can be achieved by the 
transmitter beamforming. Therefore, we conclude that trans¬ 
mitter beamforming is an optimal scheme for the MISO full- 
duplex two-way channel. In Lemma 2 we derive the closed- 
form optimal weights for transmitter beamforming. 

Lemma 2: For node i in the MISO point-to point full-duplex 
wireless network with the transmit power constraint P t and 
complex channels hu,hij,i,j £ {l,2},i j, the optimal 
beamforming weights have the following form: 


\fzi(Ci + el) ^hi 


( 12 ) 


h-ij(Ci + el)- 1 hij 

where Ci = Diag(|/i^| 2 ,..., 1 2 ), constant z t is within 

the range 0 < Zi < Pj||ft,jj|| 2 and I denotes the M x M 
identical matrix. For a fixed z, t , nonnegative constant e is 
adjusted to satisfy the transmit power constraint ||iUi|| 2 < Pi. 
Specially, e = 0 if 


PiihljCr'h^f 


Zi < 


03) 


t -2 l 
n ij ° i n i] 

Proof: The optimal beamforming weights can be obtained 


by solving problem (111 with the rank-one constraint Q ■ = 
WiW J as follows: 


min w\ C, iv, 
subject to \w\hif 2 = z it ||uti|| 2 < Pi- 


(14) 



Fig. 3. The achievable rate regions for the symmetric MISO full-duplex 
two-way channels with /3 = —40 dB, ||Iii 2 || = ||/i 2 l|| = 1, Pi = P 2 = 1. 
As plotted for comparison is the half-duplex TDMA achievable rate region. 


The above problem has the general closed-form optimal so¬ 
lution © (see details in |20j). Without the transmit power 
constraint ||«!i|| 2 < P,, problem (14» has the following 
optimal solution (shown in ]20)) 


yfziC i hi 

htC-'hi 


Combining with the condition (13 1 , we obtain 

it n~ 2 i 


112 


ZihijCi hij 




L hij) 2 


< Pi. 


Hence, we conclude that e = 0 under the condition ( fl3| l. ■ 

We remark that the optimal beamforming weights for node 
i are closely parallel to the cross-node channel hij, beam¬ 
forming the signal of interest at node j. While the transmit 
front-end noise corresponding to the stronger self-interference 
channel is largely suppressed via the matrix (Ci + el) -1 - 


V. Numerical Examples 

We present the achievable rate regions for the MISO full- 
duplex two-way channels in Fig. [3] where the channels are 
symmetric i.e., /r 12 = h 2 i and hn = h 2 2 . Each node is 
equipped with M = 3 transmit antennas and single receive 
antenna with the transmit power constraints Pi = P 2 and 
the receiver noise 0 % = 1. Define 7 = 

(in dB) to represent the ratio of the self-interference channel 
gain to the cross-node channel gain. Note that 7 can be 
reduced by the passive suppression (3j. The transmit front- 
end noise level is fixed with /3 = —40dB. Each colored line 
represents the Pareto boundary of the achievable rate region 
for the channel with corresponding 7 . We conclude from the 
numerical results that the achievable rate region shrinks as 7 
varies from 20dB to 60dB. However, the full-duplex channel 
always outperforms than the half-duplex TDMA channel if the 
optimal beamforming is employed. The extreme points A, B 
of the rate regions on the axes represent the maximum rates in 





















Fig. 4. The geometric comparison of the full-duplex optimal beamforming 
with the zero-forcing beamforming 


Fig. 5. The achievable rate regions for the symmetric MISO two-way 
full-duplex channel with = —60 dB. Circles denote the respective ZF 
beamforming rates for 7 = 20 dB ,7 = 40 dB ,7 = 60 dB and the ideal 
full-duplex network. 


the case that only one-way of the two-way channel is working. 
It follows that the points A, B are only determined by the 
transmit power constraints Pi. The ideal MISO full-duplex 
two-way channel sets the outer bound for the achievable 
rate regions of all channels, doubling that of the half-duplex 
TDMA channel. 

The rate pair with the maximum sum-rate corresponds to 
the point with Ii\ = R 2 on the Pareto boundary and can be 
achieved by certain weights in the set of optimal beamforming 
weights defined in Lemma 2. In Fig. [4] we compare the 
optimal weights OF which corresponds to the maximum sum- 
rate with the zero-forcing (ZF) beamforming weights OZ. 
For simplicity, we assume all channel vectors to be 2 x 1 
real vectors with hn = h 22 and hi 2 = h 21 . Assume the 
transmit power constraints P\ = P 2 = P- Then all possible 
beamforming weights are contained in the disc with radius 
y/P- OZ restricts the transmit signal orthogonal to the self¬ 
interference channel whereas OF is not orthogonal to h t! 
but has greater length of projection on the cross-node channel 

hij. 

In Fig. [5] we evaluate the performance of the ZF beam¬ 
forming for the same full-duplex channels as in Fig. [3] but 
with /3 = —60 dB. Comparing the channel with the same 7 
in Fig. [3] and Fig. [5] the achievable rate region is increased 
due to reduction of j3. The circles, which represent the rate 
pairs achieved by the ZF beamforming, are below the corre¬ 
sponding Pareto boundaries except for the ideal full-duplex 
channel. It follows that the ZF beamforming is not optimal 
for the MISO full-duplex channel in presence of the residual 
self-interference. As shown in Fig. [4] the ZF beamforming 
generates no interference at the unintended receiver if the 
interference equals to the projection of the transmit signal on 
the interference channel. However, for the full-duplex case in 
Q, the residual self-interference signal statistically depends on 
the transmit signal power rather than being the projection of 
the transmit signal on ha- Therefore, the ZF beamforming is 
inefficient in the suppression of the residual self-interference. 


With 7 decreasing, the residual self-interference signal is 
gradually weaker and thus the ZF beamforming is closer to 
the optimal beamforming and is exactly optimal in the ideal 
full-duplex network. 


VI. Conclusion 

We considered the MISO point-to-point full-duplex wireless 
network. We derived the achievable rate region and the charac¬ 
terization of the Pareto boundary for the MISO two-way full- 
duplex channel in presence of the transmit front-end noise. 
Using the decoupling technique and SDP reformulation, we 
proposed a new method to obtain the entire Pareto boundary 
by solving a family of convex SDP problems, rather than the 
original non-convex problems. We showed that any rate pair 
on the Pareto boundary can be achieved by the beamforming 
transmission strategy. Finally, we provided the closed-form 
solution for the optimal beamforming weights of the MISO 
full-duplex two-way channel. 


Appendix 


Proof of Theorem 1 : We prove Theorem 1 by the primal- 
dual method. Note that problem {U} is feasible and bounded. 
It follows that its dual problem is also feasible and bounded 
119). Assume Q* is an optimal solution for problem ( 111 . 
From (T9), problem (fTT| has the dual problem as follows: 


min \\Zi + A 2 P 

Ai, A2 

subject to Z = Ci — \\Ai — A 2 1 f 0. 


(15) 


where P = tr (<?*)• Assume ((A*, A 2 ), Z*) are the optimal 
solutions for We denote the rank of Q by r. We 

assume r > 1. Following that Q* is positive semi-definite, 
Q* can then be written as Q* = VV via the singular-value 
decomposition where V € C Mxr . 
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Next, we consider the following two linear equations de¬ 
fined by Ai and V : 


tr{V f A t VX) =0 

tr(V t V\X') = 0 


where the unknown matrix X £ H r contains r 2 real-valued 
unknowns, that is, ' for the real part and ' for the 
imaginary part. 

The linear system © must have a non-zero solution, 
denoted by X*, since it has r 2 unknowns where r > 2. 
By decomposing the Hermitian matrix X* , we obtain X* = 
UHXJ\ where U is an r dimensional unitary matrix and £ 
is the diagonal matrix, £ = Diag(<Ji,..., ay). Without loss of 
generality, we assume \<ji\ > |er 2 | • • • > |oy|. Non-zero matrix 
X* has at least one non-trivial eigenvalue, thus \<ti\ > 0. 
Next, we construct a new matrix as follows: 

Q* {l) = V(I - —X*)Vl (17) 


Note that I — J X“ A 0. It follows that Q*v t ’ s positive semi- 
definite. Next, we show that is also an optimal solution 
for problem ( 111 . Note that Q*^ is optimal for problem (11 1 if 
and only if (t/)' j , (A^, A 2 ), Z*) satisfies the KKT conditions, 
including the primal feasibility, the dual feasibility and the 
complementarity 1211. As ((A*, A 2 ), .Z*) is unchanged, the 


dual feasibility is automatically satisfied. Therefore, we need 
only to prove the primal feasibility and the complementarity 

°f <?(!)• 


Q(i) is a feasible solution for problem ( 111 , since the 
following two equations hold for Q*i), 


tr(4Q( 1} ) = - —X*)V') 


dl 


1 


= tr(V t AiVI) - —tr(V'AiVX*) = z i: 

o-i 


tr(Ql 1) )=tr(V(I-j-X*)V') 

= tr (vWl) - -trfytyjf*) = P. 

cr i 

To show the complementarity, note that tr (Q*Z*) = 
tr (V f Z*V) = 0 and V f Z*V t 0 implies that V f Z*V = 0. 
It follows that 


1 


tr{Q* a) Z*) = tr {V(I -X*)V T Z*) 

K ' (7 1 

= tr ((I - —X*)V^Z*V) = 0. 


Therefore, Qm is an optimal solution for problem ( 111 . 
Furthermore, the rank of QfU) is strictly smaller than r since 
rank(Q( 1 )) = rank(J — f-X*) < r. 

We can repeat this process as Q*, Q( 2 ), • • •, until 

ranl<(Q) A .j) < \/2. In other words, the rank of the optimal 
solution can be strictly decreasing to rank(Q* fc ^) < \J2, that 
is, rank(Q( fc) ) = 1 . 
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